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Abstract 

We numerically test the correspondence between the scaling limit 
of self-avoiding walks (SAW) in the plane and Schramm-Loewner evo- 
lution (SLE) with K = 8/3. We introduce a discrete-time process 
approximating SLE in the exterior of the unit disc and compare the 
distribution functions for an internal point in the SAW and a point at 
a fixed fractal variation on the SLE, finding good agreement. This pro- 
vides numerical evidence in favor of a conjecture by Lawler, Schramm 
and Werner. The algorithm turns out to be an efficient way of com- 
puting the position of an internal point in the SAW. 



1 Introduction 

Schramm-Loewner evolution (SLE) is a family of random processes on confor- 
mal maps, which gives rise to a one-parameter family of measures on curves 
satisfying conformal invariance. It was first introduced in pQ to describe the 
scaling limit of loop-erased random walks, but it was soon found to corre- 
spond to a large family of geometrical objects defined in the context of lattice 
models. The latter include interfaces in critical models — e.g. the interface 
between phases in the Ising model or the boundary of the percolating cluster 
in critical percolation — and walk models, such as the self-avoiding walk 
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in the scaling limit. A wealth of results have been obtained for SLE in the 
chordal geometry, which involves curves starting and ending on points ly- 
ing on the boundary of some connected domain, and in the radial geometry, 
where one of the points is on the boundary and the other is in the bulk. 
Much less has been proved or checked in the whole-plane geometry, where 
both the starting and ending points lie in the bulk, such as the points and 
oo in C. 

In this paper we consider SLE in the whole plane and compare it with 
the planar self-avoiding walks (SAW). The correspondence between the two 
models has been conjectured by Lawler, Schramm and Werner in [2], on the 
basis of restriction covariance and conformal invariance. In particular, we will 
compare the scaling forms of the distribution functions in the two models. 
For this purpose, one has to be careful about which points along the curves 
have to be considered. On the SLE side, a suitable parametrization must be 
chosen; we will rely on parametrization by fractal variation, and will pick the 
point at a fixed value of this parametrization. On the SAW side, we are going 
to consider an internal point deep inside the walk, so as to avoid finite-chain 
effects. Thanks to the scaling relation between the fractal variation and the 
natural parametrization of SLE, and to the fact that the discrete chains are 
sampled independently, this turns out to be a potentially efficient method for 
computing the distribution function of an internal point in the SAW, since 
it produces independent samples of length L in time 0{L^'^) with room for 
easy improvement. 

The plan of the paper is as follows. In section [2] we quickly review the 
relevant facts and definitions about Schramm-Loewner evolutions and self- 
avoiding walks; in section [3] we introduce and discuss the discrete process we 
will be simulating in the whole plane; in section H] we define the distribution 
function and its conjectured scaling behavior; in section El we present our 
numerical results; the appendix deals with the scaling of the average step 
length and the hull size in discrete SLE. 

2 Definitions and background 

SLE is a widely studied and reviewed subject, see e.g. [3], H] for an intro- 
duction. Here we will only recall the main ingredients, in order to make this 
paper as self-contained as possible. 

SLE is a stochastic differential equation describing the evolution of a 
parametrized family of conformal maps gt in a domain T>. The actual form 
of the equation depends on the geometry, but in general it is written in terms 
of a function taking values on the boundary of V, called driving function. The 
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driving function is rescaled Brownian motion y/K,Bt living on the boundary; 
K is a positive constant. The maps gt map some (t-dependent) subset of V 
onto so that one can define the growing hull Kt as the domain of gt, i.e. 
the set of points for which the differential equation still has a solution up 
to time t. The hull turns out to be generated by a curve ■jit), i.e. it is the 
union of the image of the curve with the interior of any loop it has closed. 
The curve — also called the trace of the process — is the pre-image of the 
driving function under gf. 



As the variance k of the Brownian driving function varies, the properties 
of 7(t) change dramatically, and very different models have been proved 
or conjectured to be described by SLE for some value of k in the scaling 
limit. SLE has proved a useful tool for addressing very diverse problems and 
questions (see for instance [51 E] and references therein). 

Interesting issues arise when one considers the reparametrization of 7(t) 
[7]. The SLE equation generates the hulls with a natural parametrization, 
which corresponds to a linearly growing capacity (which is a property related 
to the expansion of gt around some special point and is a measure of the 
"conformal size" of the hull, see ^ for a precise definition). As long as one 
is interested in parametrization-independent observables this is not an issue. 
But we will be considering the distribution function of the point 7(t) for some 
fixed t, so the choice of parametrization is clearly essential here. One way of 
dealing with this problem in the context of a numerical work is to generate the 
curves already with the correct parametrization — here correct means the one 
corresponding to the parametrization by length of the SAW. The technical 
problem of producing discrete SLE traces parametrized by length will be 
treated elsewhere. In this paper instead we will take advantage of a technique 
introduced by Kennedy [9], which consists in defining a suitable notion of 
length along the discretized curve, and using it as the parametrization (see 
section [H]) . 

A self-avoiding walk u on the square lattice with a fixed number of steps 
is an ordered collection u = {ujq, . . . ,ujn} with a;, G Z^, such that 
\uji — cjj-il = 1 and a;, 7^ ujj for i ^ j. The SAW model is the uniform mea- 
sure on these objects, i.e. on all non-intersecting A^-step nearest-neighbor 
walks 

A numerical check of the equivalence between the distribution functions 
of SLE and SAW in the half plane has been carried out in [9] , with positive 
results. 



gt (7(t)) = v^Bt. 
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3 Discrete whole-plane SLE 



Simulation of SLE in the half plane EI is usually based on the discrete process 
introduced in [T2], where it is also proved that it converges weakly to SLE. 
Here we introduce a whole-plane version of this discrete process. 

The idea (see also [TUj) is to write the Loewner map as a composition of 
(finitely many) conformal maps, chosen from a simple parametrized family. 
Time is partitioned by letting = to < < • • •• The driving function 
is defined in such a way as to be equal to y/^Bt^. at the special times tk 
and constant in between. This amounts to approximating the Brownian 
motion with a piecewise constant function, so that the incremental maps 
Gk = 9tk ° 9r^_ obtained by solving the Loewner equation with constant 
driving function up to time tk — tfc-i. 

Let us denote = \fi^Bt. Consider the radial SLE equation 



which describes a hull starting from exp (i^o) and growing towards the origin 
inside the unit disc D. In analogy with the half-plane case, the hull Kt at 
time t is defined as the set of those points in © for which the differential 
equation does not admit a solution that exists up to time t. By composition 
with the complex inversion {z ^ l/z) one can define the map 



which now describes a hull growing in C \ D from exp (— i^o) to oo. The 
new map in ([3]) happens to satisfy the same equation ([2]), again with initial 
condition gQ{z) = z, but now with driving function —^t, which has the same 
law as if started from = 0. This process is called radial SLE growing to 
infinity. Notice that gt does not describe a curve truly in the whole plane, 
since the unit disc is a forbidden region. However, we are going to focus on 
the large-scale regime in the following, and expect the cutoff at length scale 
1 to be irrelevant in this limit. 

The discrete process is based on the discretization of the evolution of the 
inverse map g^^. This map grows the hull at time t, i.e. it maps C \ © onto 
the complement of the hull in C \ D. The atomic step of the discretization 
will be performed by a map — which we are going to call incremental — that 
grows a slit (a small radial segment) out of the unit disc. Obtaining this map 
is simply a matter of solving ([2]) in the special case when the driving function 
is constant, thus finding the Loewner map corresponding to a straight line 







(3) 
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growing towards infinity; the incremental map is its inverse. The resuh i^ 



"^'^'^ = 2rt [^^ + ~ - + 1) ((^ + 1)' - 4e-*^)'/'] . (4) 

The greater the time t is, the longer will be the line grown by (pf In the lan- 
guage of potential theory, one says that the hull grown by (pt has logarithmic 
capacity t. 

Brownian motion driving the evolution lives on the boundary of the stan- 
dard domain, which is the upper half plane in the chordal case, the unit disc 
in the radial case and the complement of D in C in the case at hand. Notice 
that, by translation invariance, chordal SLE in H growing from to oo is 
chordal SLE growing from to cxd translated by ^o, so that in simulating the 
discrete process one usually translates back to the origin after each itera- 
tion. Translations are replaced by rotations in the radial case. Therefore, at 
each step we will want to rotate exp (i^J back to 1. The whole discretized 
hull will be generated by alternately composing an incremental map with a 
rotation. Let us call the logarithmic capacity of the incremental map at 
the k-th step (i.e. the time at which the k-th map is evaluated), and corre- 
spondingly let 6k denote the angle of the A;-th rotation. The incremental map 
and the rotation themselves will then be denoted and i?^^. respectively, 
where Rs^i^z) = zexp{i6k). Let 7^ be the image of 1 under the composed 
map at step k 

7fc = -R^i o o ■ • ■ o o 0Afc(l). (5) 

Notice that the order in which the maps are composed is the opposite as the 
usual one; this is essentially a consequence of the fact that we are discretizing 
the inverse Loewner map. The approximate curve we are interested in is 
embodied by the collection of points {jk}- 

Refer to figure [H The first map grows a slit based at point 1, which is 
then rotated to exp(i5fc). The second map again grows a slit based at 1, 
so that the old slit will be sent away from the unit disc, and will be based 
somewhere between 1 and the tip of the new slit (notice that the old slit will 
not in general retain its rectilinear shape). Then the whole hull is rotated 
by 6k-i and another slit is grown. When the last composition with Rs-^ o 
is reached, the k-th point on the trace is found. The process is repeated for 
each point 7^ that has to be computed. 



^ Some care must be taken in choosing the right sign before the square root and the 
position of the branch cut of the square root itself, so that a point outside D be mapped 
to a point outside D. The conformal map in ^ aheady appeared in the hterature about 
diffusion hmited aggregates , where it was used to represent the attachment of a single 
grain onto the cluster. 
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Figure 1: (Color online) A visualization of the discrete process. The leftmost 
arrow corresponds to the action of i?^^, o 0^^, the second corresponds to 
-R^fc-i ° <^Afc_i, and so on. The blue cross is the point 1, where the next slit 
will be based. 



The resulting trace depends on the choice of the time-like parameters 
{Afc} and the space-like parameters {Sk}- In analogy with the chordal half- 
plane case we choose to draw each Sk as a Bernoulli variable in the set 
{ — -v/^Afc, -f-A/KAfc}. This choice amounts to approximating the Brownian 
motion with a piecewise constant function, the relation between and 
6k reproducing the well-known space-time scaling of Brownian motion with 
variance k. One has additional freedom in choosing the time intervals {A^}. 
Changing the latter corresponds to a reparametrization of the resulting trace. 
The uniform partition = A yields the parametrization by capacity, i.e. the 
natural SLE parametrization where the hull increases its capacity linearly in 
time. The choice of parametrization is crucial for our purposes, since we are 
interested in the spatial distribution function, which is a parametrization- 
dependent observable. In order to reproduce the correct parametrization 
(that corresponding to the natural parametrization of the supposed scaling 
limit of self-avoiding walks) we use the following method. First, we fix a 
scale A. Then, at each step, we compute the fractal variation of the curve as 
follows. Let ko = 0. Times {ki}-^^ ^ are defined recursively: given fcj, fcj+i 
is the first time after ki such that l^ki+i ~ TfeJ > A. The fractal variation at 
step k is defined as 

varA({7o,...,7fc}) =nA'^^. (6) 

where n is the largest integer such that kn < k, and df is the fractal dimension 
of continuum SLE function of k [14J 

df = l + l. (7) 

The growth process will be stopped when var;^ reaches a fixed value T. This 
procedure has the advantage of being independent of the original param- 
etrization of 7. Nonetheless, it is sensitive to discretization problems, for 
instance when the steps taken by the trace 7^ become too wide with respect 
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to A. A uniform partition At = A causes the average step length at step k, 



h = (|7fc ~ lk-i\), to diverg^ as the number of steps k grows (the average 
is over all realizations of the process up to step k). To avoid approximation 
problems we choose the time intervals A^ in a non-uniform fashion, having 
them scale as 

A, ~ k-' (8) 

so as to compensate for the divergence of the average step length (see the 
appendix for details). We can then define / = 4, at least for k large enough. 
The discrete process is defined operatively as follows: 

Fix the three scales /, A and T; start with k = 1 
Choose Afc according to (IHl) 
Draw 6k uniformly in | — -y/ztAfc, +-\//cAfe} 
Calculate 7^ as in 
Measure the fractal variation 

If varA ({70, • • • , 7fc}) > T stop, otherwise increase k and return to 
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Notice that the time when the fractal variation reaches T depends on 
the scale A it is measured at, which should be sent to in order to obtain 
the true fractal variation. Of course, the curve obtained by means of the 
discrete process defined above is not really fractal at all, but it displays fractal 
properties only at a large enough scale. Therefore, values of A much less than 
the average step length are expected to yield the trivial parametrization. 
On the other hand when A becomes comparable to T^^'^f , vaix suffers from 
rounding problems. We will keep A between these two cutoffs, and study how 
results depend on this choice. 

An example of the composed map after a few iteration is in figure El 



4 Distribution function 

Let us consider the M-th point in an A^-step SAW u. The probability distri- 
bution PN,M{f) of this internal point (the probability that it lies at site r on 
the lattice) has the following scaling form in the limit N, M ^ 00, r — > 00, 
with M/N fixed: 

PNMi.f) ~ ^/SAW (^P, ^ (9) 

^ This is true in the whole plane. The opposite happens in the half plane, where the 
average step length converges to zero. 
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Figure 2: A visualization of the composed map in ([5]) after 4 iterations. The 
curves are the images of circles and radii outside the unit disc. 

where p = 1?"/^^/ 1 and = {\u}m — ^of) is the correlation length after M 
steps. The average (■) is on the self-avoiding walk ensemble. The univer- 
sal function /saw is the renormalized distribution function we are going to 
compare with its SLE analogue. It depends on the ratio M/N — i.e. on 
how much the M-th point feels the finiteness of the chain. Since we want 
the distribution of a point a finite distance away from the origin in a truly 
infinite curve, we should take M/N — > 0. 

Let us then define the corresponding quantities for the SLE approximated 
trace {7^} {k = 0, . . . , Ny), where A^x is the number of points computed on 
the trace up to the stopping time when the fractal variation has reached T. 
We are not considering the full distribution P{z), since the discrete process 
as defined above explicitly breaks rotational invariance — it has 1 as a special 
point. Instead, consider the probability density Prir) that the point 7^ has 
modulus r when k = N^, then we will suppose it has the following scaling 
behavior in the long-chain limit, i.e. when T ^ 00 and r ^ 00 with step 
length / fixed: 

Pr{r) ~ ^ /sle(p) (10) 

ST 

where p = r/C,r, and is the correlation length 

er = {\lNrf) (11) 

(the average (■) here is over all realizations of the driving function). Implicit 
in this conjectured behavior is the assumption that the presence of the unit 
disc as a forbidden region be irrelevant in the long-chain limit defined above. 
Scaling form (JTOl) for the distribution function will be verified a posteriori. 
Notice that of course /sle implicitly depends on k. 
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In order to quantitatively compare the distributions we will focus on 
their moments. It is convenient to introduce an infrared cutoff Pmax — i-e. a 
window for computing the moments — since deviations from scaling are more 
pronounced in the large-p regime. Suppose we have sampled TV instances Zi 
= 0, . . . , A/" — 1) of one point on the chain — be it an internal point inside 
a SAW or the point at a fixed fractal variation on the approximated SLE 
trace. We will compute the following quantities: 

M^k = ^ , (12) 
[E \M \ 

where the sums Yl' restricted to the window 

< Pmax- (13) 



\Zi\ 



We will choose pmax = 3 in the following, which only leaves out the tail of 
the distribution (see figure HI). 



5 Numerical results 

We simulated an ensemble of 10^ self-avoiding walks of length = 100 000 
using the pivot algorithm [THl IS]. We considered an internal point with 
M = 8000 (detailed analysis of the systematic error due to the finiteness of 
M/N = 0.08 shows that it is negligible when compared to the deviations 
in discrete-SLE data due to the finiteness of A and T). The discrete SLE 
process was simulated for several values of A and T, generating ~ 100 000 
independent samples for each choice. In the following we are going to fix the 
average step size / and measure everything else in units of /. 

The time needed to compute the fc-th point along the trace through ([5]) 
is proportional to k. Therefore, generating a A^-step chain requires a time 
of order A^^, which can be interpreted by saying that the time-per-point is 
O(A^). As shown in the appendix, partitioning time as in ([H]) causes the 
fractal variation to scale as T ~ N'^f . All together, the algorithm described 
here generates chains of length T in time 0{T'^^'^f), which is 0(T^/^) for n = 
|. The time-per-point can be probably further improved by approximating 
the incremental map by its truncated Laurent series (see [10], where the 
numerical analysis is carried out in the half-plane geometry, suggesting that 
the time-per-point is 0(A^°''^)). We shall not do that here, but this suggests 
that this algorithm might get close to generating independent samples of 
length T in time O (T), and even faster for higher values of k,. Producing 
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Figure 3: (Color online) The number of segments of length A needed to cover 
the trace (red crosses), for k = | (left) and k = 6 (right). Dashed green lines 
represent the theoretical slope corresponding to fractal dimension 1 + k/8. 
(Error bars are smaller than the crosses representing data). 



100 000 samples with T = 400 took about 1000 hours on an Intel Pentium 4 
with 1.80 GHz CPlH. 

As a preliminary test, we checked that the discrete SLE approximants 
approached the expected fractal dimension (JTj), since knowing its precise 
value is crucial whe; computing the fractal varil^iorfl as in ([6]) . For increasing 
values of A we measured the number n{X) of segments of length A that are 
needed to cover up the entire trace, in the same fashion as when computing 
the fractal variation. The expected behavior is n{X) ~ X~'^f . We checked this 
for several values of k, finding good agreement. Figure [3] shows the results 
for two different values of k. For small A one sees the crossover to the true 
fractal dimension of the discrete trace, which is of course 1. 

We are now ready to check whether /saw = /sle for k = |. Figure H] is a 
comparative plot of the renormalized distribution functions for the SAW and 
the discrete whole-plane SLE, for two values of the variation, T = 400, 800 
(an averaging procedure has been adopted here, in order to smoothen oscil- 
lations due to the lattice). The distributions all fall onto the same universal 
curve, apart from corrections to scaling in the extreme regimes. To quanti- 
tatively check that this is the case, the moments ffT^ have been computed 
for each value of T and for several values of A; they are reported in table [T], 
together with their SAW values. Accordance for low-order moments is more 

^ This is the performance of a non-optimized code. 

^ Notice that the fractal dimension does not depend on the parametrization. Actually, 
we reverted to parametrization by capacity when testing df — i.e. we stopped the discrete 
process at a fixed number of iterations, without computing the variation. 
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Figure 4: (Color online) The renormalized distribution functions of SLE with 
T = 400 (red squares) and T = 800 (green circles), and that of a point inside 
a SAW (blue triangles). The two plots display the same data; the one below 
has logarithmic scale on the y-axis. 



easily established, while inspection of high-order ones helps in recognizing 
systematic deviations. Let us fix T = 400 first. Values for A = 40 suffer 
from severe deviations due to the coarse grained nature of the procedure 
used to compute the fractal variation. On the other hand, a systematic drift 
in A is present for smaller values. The best compromise seems to lie in the 
middle; we will then fix A = 20 and increase T. The moments for T = 800 
all lie within two standard deviations from the SAW values, apart from the 
highest-order one which is close. Deviations are apparent, the SAW values 
being systematically larger than the others, but again SLE moments keep in- 
creasing as T is doubled, thus approaching the expected values. As expected, 
these systematic corrections get larger if the cutoff Pmax is increased. Con- 
vergence for p = 3.5 (see M12 in the bottommost line in table [T]) is slower 
but still consistent. 
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T = 400 T = 800 

X^4 A = 10 A = 20 A = 40 A = 20 SAW 



M4 1.304(10) 1.315(7) 1.320(11) 1.318(7) 1.322(9) 1.330(2) 

Me 1.970(18) 2.009(14) 2.022(25) 2.008(14) 2.030(20) 2.059(5) 

Mg 3.267(38) 3.374(32) 3.404(55) 3.352(31) 3.428(42) 3.508(10) 

Mio 5.787(84) 6.045(70) 6.11(12) 5.955(68) 6.177(93) 6.379(23) 

M12 10.74(18) 11.35(16) 11.48(27) 11.07(15) 11.66(21) 12.155(52) 

M12 15.77(32) 16.94(27) 16.79(46) 14.92(23) 17.43(37) 18.763(94) 

Table 1: The first five non-trivial moments of the distribution function for 
whole-plane SLE and SAW. T is the fractal variation computed at scale A. 
(Here Pmax = 3, except for the bottom line where Pmax = 3.5). 

6 Conclusions 

We have introduced a discrete process approximating radial Schramm-Loew- 
ner evolution growing to infinity, by the iteration of conformal maps defined 
outside the unit disc. On one hand, we have considered the distribution 
function of a point on the trace at a fixed value of the fractal variation. 
As the fractal variation reaches infinity the distribution forgets about the 
presence of the forbidden region D. On the other hand, we measured the 
distribution function of a point deep inside a whole-plane self-avoiding walk. 
When K = |, the two universal functions match, thus providing evidence 
that the scaling limit of SAW is whole-plane SLE. Moreover, computing the 
position of an internal point in a SAW by exactly sampling the discrete 
SLE process seems to be an efficient algorithm, which is open to further 
improvement. 
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Appendix 

We derive here the scaling form ([8]) . Let Ik denote the length of the A;-th step 
in the discretized growth, that is 

L = bn - 7n-l| • (14) 
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The exponential of the logarithmic capacity is a measure of the linear size of 
the growing hull. This is essentially a consequence of the Koebe 1/4 theorem, 
which bounds the size of the image of the unit disc D under a conformal map 
g in terms of 5''(0). Then, by approximating by 



In ~ |7n| - l7n-l| 

as if it grew radially, one has 

vfc = l 

(e^"-l) expfJ^A, 



(15) 



~ exp 



exp 



'n-l 



n-1 
.fc=l 



(16) 



Setting 

in f|T6|) gives 



A 



exp — - 1 ) exp {AHn-i] 
n 



where the ifn's are the harmonic numbers, whose expansion in n is 



Hr, 



In n + 7 + O I — 



n 



(17) 
(18) 

(19) 



(7 is the Euler-Mascheroni constant). By substituting in ffT^ and expanding 
both exponentials in n one gets 



Ir. 



A 



n 



+ 



[n 



1 + 



(20) 



By choosing A = 1 one finally obtains 



1 + 



(21) 



which shows that scaling A^ as in (JHl) provides an asymptotically constant 
step length. Contrary to the half-plane case — where a similar computation 
shows that an approximately constant /„ is obtained by choosing A^ = kA 
— there is no freedom left here to choose the average step length, since A has 
to be fixed to 1 in order to have the correct scaling behavior. Operatively, we 
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keep a constant for a few steps, until the step length has approximately 
reached the desired value I, and we set = 1/k from then on. 

In order to obtain a scaling relation for the average number of steps n 
needed to reach fractal variation T„, again we use the exponential of the 
logarithmic capacity as a measure of the chain size, which by ([6]) is of order 
tI/'^^ , so that for n large we have 

n j_ 

expJ]Afc~T^^ (22) 

k=l 

which, by substituting Ak given by ( fTTl) with A = 1 and again using ( fT9l) . 
yields 

n~T^. (23) 
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